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We examine the leptophilic two Higgs doublet model with fermionic dark matter, considering the 
range of experimental constraints on the Higgs sector. The measurements of the 125 GeV Higgs 
from the LHC Run-I allow us to focus on those remaining processes that may play an important 
role at colliders. We find that the leptophilic model allows for a much lighter Higgs than in other 
two-Higgs models, although discovery at the LHC will be difficult. Adding a dark matter sector 
motivated by supersymmetric extensions of the leptophilic model, we find the existing parameter 
space can accommodate constraints from direct detection and the invisible widths of the Higgs and 
Z, while also fitting the Galactic Center gamma ray excess reported by analyses of Fermi-LAT data. 
We also discuss the status of the fully supersymmetric version of such models, which include four 
Higgs doublets and a natural dark matter candidate. 


I. INTRODUCTION 

The minimal particle content necessary to realize the Higgs mechanism in the Standard Model is one scalar SU (2 )l 
doublet with hypercharge ±1/2. The discovery of the 125 GeV Higgs by ATLAS and CMS [U |2] is an important 
step on the path to understanding electroweak symmetry breaking (EWSB) and appears consistent with the minimal 
Higgs sector. However, we should keep in mind that the full scalar sector remains unknown. Therefore, searches for an 
extended Higgs sector will be an important part of the LHC Run-H. One of the simplest ways to enlarge this sector is 
to introduce additional scalar doublets, with the minimal extension being a two-Higgs-doublet model (2HDM). While 
the best known realization of these additional doublet models occurs in supersymmetry, which requires a minimum 
of two Higgs doublets and is perhaps the most popular candidate for physics beyond the Standard Model, 2HDMs 
can be found in a variety of other Standard Model extensions. For example, models of baryogenesis with two Higgses 
[SHE] can incorporate enough CP-violation to generate a sufficient baryon asymmetry, while the Standard Model 
alone cannot [T5] . 

Whereas the Standard Model Higgs has one physical degree of freedom after EWSB, in the form of a charge-parity 
(CP) even scalar, 2HDMs have hve degrees of freedom: two CP-even neutral scalars {h and H), one CP-odd neutral 
pseudoscalar (a), and a charged scalar (H^). With two doublets the fermions can couple to the scalar sector in a 
variety of ways. In this paper we examine a 2HDM in which one Higgs doublet has Yukawa couplings with both the 
up- and down-type quarks, while the other doublet couples to the leptons. These are known as leptophilic two-Higgs 
doublet models (L2HDM).^ 

Though studies of leptophilic models (including connections to dark matter) have been carried out previously |14F 
[20] , after discovery of the 125 GeV Higgs and the early measurements of its couplings to gauge bosons and fermions 
(see Refs. for experimental results and Refs. [251 - HD] for theory analysis) we can paint a much clearer picture 

of the available parameter space in 2HDMs, leptophilic and otherwise. In particular, the CP-even Higgs boson at 
125 GeV appears to be close to the alignment limit, resulting in properties very similar to those found in a minimal 
single-Higgs sector. As a consequence, processes which would be potentially significant for production of additional 
physical Higgs scalars at the LHC in a generic 2HDM are now expected to be highly suppressed. 

With this in mind, we consider the available parameter space for a L2HDM, reinterpreting Run-I searches and 
considering the prospects of future searches at the 13 TeV LHC. We further consider the possibility of coupling 
fermionic dark matter to the Standard Model through the leptophilic Higgs. So-called Higgs Portal couplings to dark 
matter are theoretically attractive [nmg and link two of the most pressing open questions in theoretical physics. 
As we shall show, leptophilic Higgs models greatly relax the existing constraints as compared to other 2HDMs, and 
suggest specific channels for further searches at the LHC. As fermionic dark matter coupled to the Higgs sector is 
realized in supersymmetric extensions, in addition to the “simple” L2HDM, we further consider a supersymmetric 
extension, which contains minimally four Higgs doublets. 

In addition to the experimental constraints on the Higgs sector and new particles coupling to gauge bosons from 
the LHC and LEP-H, we consider the leptophilic Higgs contribution of the new scalars to the {g — 2) value of the 


^ In the literature the leptophilic 2HDM is also commonly referred to as the Type-X 2HDM. In other works it has also been called “Model 
I”, “Model IIA”, and “Model IV.” 
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muon, on dark matter from direct detection experiments and from indirect detection searches. In particular, we 
discuss the dark matter interpretation of the unexplained excess of gamma rays coming from the Galactic Center, 
found in the Fermz-LAT data |43H46) . Fermionic dark matter annihilation mediated by scalars is p-wave suppressed, 
i.e. the cross section is w^-dependent at leading order. In this case, the velocity of dark matter today is far too small 
to account for the necessary cross section needed to explain the gamma ray excess. However, annihilation modes 
through a pseudoscalar proceed via the velocity-independent s-wave channel, and so has potential to fit the anomaly 
(as has been pointed out in the leptophilic context in Ref. m)- We find parameter space available to account for the 
observed excess, compatible with the LHC Higgs results, in both the L2HDM and the full supersymmetric four-Higgs 
model. The solution requires a light pseudoscalar and large tan /3, which is also the region of interest for solutions to 
the {g — 2)^ anomaly. 

This paper is organized as follows. Section [IT] gives an overview of the leptophilic two-Higgs doublet model. The 
theory along with constraints from the Higgs sector are described in detail. Section [m] introduces the dark matter 
sector of the L2HDM. The ingredients are presented along with notation and couplings. We follow with discussions of 
implications for direct/indirect detection experiments, as well as for the thermal relic abundance. Section IV discusses 
the full supersymmetric model and the key differences between it and a generic L2HDM. 


II. THE LEPTOPHILIC TWO-HIGGS DOUBLET MODEL 

All models with two or more Higgs doublets must assume some symmetry structure in the Higgs-Standard Model 
fermion couplings in order to prevent large flavor violation at tree-level. The L2HDM [TlU^n] extends the Standard 
Model to contain a Higgs sector with two complex scalars, and assumes an additional Z 2 symmetry which selects one 
doublet ($/,) to have Yukawa couplings with the leptons only, and the other ($q) to have Yukawa interactions with 
both up- and down-type quarks. We first review the notation of the L2HDM, and then consider the experimental 
constraints on this model. 


A. Higgs Sector 


Our two SU(2 )l doublet scalars we assume have equal and opposite hypercharge, being the doublet with —1/2 
and with -1-1/2. Allowing for the possibility of charge-neutral vacuum expectation values (vevs) vq and we 
can expand the fields in terms of the neutral CP-even, neutral CP-odd, and charged components: 


4>q 



(vQ + hQ+ lag) ) ’ 


= 




hr — 



( 1 ) 


Assuming both that CP is conserved in the Higgs sector and that CP is not spontaneously broken when the electroweak 
symmetry is broken, and imposing a softly-broken Z 2 symmetry under which is odd and 4 >q is even, the most 
general scalar potential is: 

+ y + y I^qI^ + + A4|$q$lP + |y + h.c.j (2) 

where all coefficients are real. Experimentally viable EWSB occurs for the region of parameter space where this 
potential is minimized by 

(^q) = ^ ^ ^ > (®i) ^ 

Requiring stability of the vacuum allows us to solve for the mass parameters fiq and g,L- We define the angle /3 by 
the ratio of the vevs tan(,S) = vqIvl, and the vevs themselves satisfy Vq + v‘1 = v‘^ = (246 GeV)^. We will typically 
be interested in the parameter space where tan/3 > 1, which will result in the physical Higgs fields having increased 
couplings to leptons and suppressed couplings to quarks. 

After EWSB, three fields are eaten by the W and Z gauge bosons, leaving five physical degrees of freedom: two 
neutral CP-even states, a neutral CP-odd state, and a charged state. The angle a rotates the CP-even states into 
the physical mass eigenstates H and h, while the angle /3 rotates the CP-odd state and the charged states into the 
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TABLE I: A summary of couplings between the Higgs scalars and Standard Model fermion and gauge boson pairs. 
There is an implied ^ multiplying the fermionic table entries, while the gauge boson entries should be multiplied 
by their Standard Model couplings. The pseudoscalar a has a + cot f3 prefactor in the coupling to u-type quarks and 
a — cot P to d-type. The right table shows the couplings in the alignment limit cos (/3 — a) = 0. 


physical states a and H^: 


a = — sin/3 Oi + cos/3 ag 

= — sin /3 + cos /3 hg (4) 

h = — sin a/iL + cos a/iQ 
H = cos a hz, + sin a/iQ, 


After diagonalizing the mass terms in the potential, we find the physical masses take the form: 




rnjj± 


^ (AlC^ + AaS^) ± - ^^^^34552/?) 


- 


S2P 

^ - y (A4 + As) = y (A4 - As), 

S2^ 2 2 


( 5 ) 


where we employ the notation Sx = sin a;, Cx = cosx, and tx = tanx, and A 34 S = A 3 + A 4 + As. Note that it is 
possible to choose these four masses independently in a general 2HDM, though the EWSB precision measurements 
do require a relatively small mass splitting between the iJ and (see e.g. Ref. [48]). This differs from the minimal 
supersymmetric Standard Model (MSSM), where the quartic interactions are set by gauge-couplings via T>-terms, 
with 


2 2 2 2 2 

^MSSM _ \MSSM _ 9l + 92 ^MSSM _ 92 ~ 9l ^MSSM _ ^MSSM _ q / 0 '| 

1 '2 ^ 5 '3 ^ ^ ^ 2 ^ 

Here gi and are the hypercharge and SU{2)l gauge couplings, respectively. A consequence of this is that the 
masses of the a, and either the h or the H are close in mass (at tree level). In particular, it is difficult to separate 
the mass of the pseudoscalar and one of the CP-even scalars in the MSSM at tree-level. 

Assigning the left-handed lepton doublet and the $z to be odd under some Z 2 , the Yukawa interactions are 
constrained to be (suppressing flavor indices) 


>CYuk 2 —VeLL^LeR — VuQl^qur — ydQL^qdR + h.c. 
After EWSB, the couplings of the physical Higgs fields to the fermions is given by 
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( 7 ) 

( 8 ) 


The dependence of the couplings between the Higgs fields and the Standard Model fields relative to the single-Higgs 
value is summarized in Table |lj As tan /3 is increased, the couplings to the leptons for the non-Standard Model-like 
scalars are increased, while the quarks are decoupled; thus in this limit the additional fields are truly leptophilic. 

In the case of a generic 2HDM, there is no reason a priori to prefer that either the heavier H or lighter h is the 
Standard Model-like Higgs discovered at 125 GeV. However, in the MSSM, the masses of rriRi, mn, and rria are 
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constrained to be of the same order (modulo splittings from EWSB), and thus assigning H as the 125 GeV Higgs 
would imply relatively light charged Higgses, which would result in large flavor-violating decays in the quark sector 
[45] , Thus is it typical to identify the lighter h as the 125 GeV Higgs, and place the H, a, and heavier. As 
we will show, one of the interesting features of the L2HDM is the relaxation of these bounds, widening the available 
parameter space to the point where the 125 GeV Higgs can be identified as the heavier CP-even Higgs. 

The measurements of the 125 GeV Higgs at the LHC implies that the couplings of the CP-even Higgs to weak 
bosons must be very close to the single-Higgs Standard Model values. Similarly, couplings to fermions, though less 
well constrained, are also consistent with the Standard Model prediction, though this is primarily driven by the 
contribution of the top-Higgs coupling to the loop-induced and hgg widths rather than from direct measurement 
of associated production. In Figure [T] we show the best-fit contours to the LHC Run-I data as a function of the 
parameters tan/3 and cos/3 — a for the L2HDM, assuming tree-level couplings.^ For notation simplicity however, we 
will in general refer to the 125 GeV Higgs as the h. The experimental data indicates that the 125 GeV h is close 
to the alignment limit, where cos(/3 — a) ~ 0. Note that had we assigned the H as the 125 GeV particle the only 
replacement in this fit would have been cos(/3 — a) O sin(/3 — a). In this work, we will assume that the h lives exactly 
in the alignment limit; though minor relaxations of this assumption within the experimental best-fit regions will not 
significantly change our conclusions. 

In the decoupling limit where rria, , and mn are large rnz), the couplings of the lighter Higgs automatically 
become Standard Model-like. Alignment is possible without decoupling (see e.g. Ref. [37]), but one may need to appeal 
to some mechanism in order to explain the observed similarity of the 125 GeV Higgs to the Standard Model predictions. 

The 125 GeV h can decay into a pair of pseudoscalars, if the a is kinematically accessible {i.e. TOq < 62.5 GeV). In 
the alignment limit, the relevant coupling is 

ghaa = g (^1 + A 2 — 2 A 345 ) -I- Ai -I- A 2 -f 6 A 345 — I 6 A 5 ) , (9) 


where A 345 = A 3 -I- A 4 -I- A 5 . This will result in a —>■ aa partial width of 
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( 10 ) 


Given the experimental constraint on the total width of the 125 GeV Higgs T/j < 17 MeV [SD], and the branching 
ratio into non-Standard Model channels must be < 50% [50H54] this coupling is limited to be \ghaa\ ^ (H GeV) x 
(1 — However, if we instead assume that this is the only significant contribution to the Higgs width on 

top of the Standard Model decays, then new physics contribution to the width can be no larger than ^ 5 MeV, and 
the limit becomes \ghaa\ ^ (8.0 GeV) x (1 — Assuming supersymmetric-like couplings, as in Eq. ^ 

results in ghaa = (—34 GeVcos^ 2/3) x (1 — which is well above the minimum for tan/3 not near 1. 

As a result, without significant loop-corrections to the couplings, supersymmetric-like models are constrained to have 
ma > 62.5 GeV. 


B. Charged Higgs Constraints 

At this time, the best direct limits on the charged Higgs scalars come from LEP searches |SS], which place a lower 
limit of mhi± > 92 GeV. The majority of LHC searches for charged Higgses rely on the coupling to top quarks 
[56|-|59|, which is tan/3 suppressed in the L2HDM.^ Direct charged Higgs pair production followed by decays into 
T+r“ -I- is searched for, but existing limits are not yet sensitive to new scalars with masses at the LEP-H bound 
m- As a result, no limits from the LHC can be set on charged Higgs from the L2HDM for tan /3 > 1 as of yet. 

Indirect limits can be derived from precision measurements of rare decays such as B ^ 57. In a Type-H 2HDM 
(as found in the MSSM) where one Higgs doublet couples to up-type quarks and the other to down-type quarks and 
leptons, the physical charged Higgs couples to up-type quarks as cot/3 and to the down-type quarks as tan/3. In 
contrast, the couplings of the to both up and down-type quarks in the L2HDM goes as cot/3 (see Table and 
Eq. §). 


^ Fit technique is motivated by Ref. m- 

^ The exception is Ref. [60], which uses — WZ couplings, but the cross section to which this search is sensitive is too large to place 
significant constraints on the 2HDMs considered in this paper. 
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cos{f3—a) 

FIG. 1: 1(7 (green) and 2cr (yellow) confidence fit of LHC experimental results pTI - B5| to the L2HDM in the 
cos(/? — a) — tan/i plane, assuming tree-level couplings and uncorrelated errors in the experimental measurements. 


Therefore, in the Type II/MSSM case the H^ud coupling cannot be fully suppressed in any limit of tan /3, forcing the 
to be heavier than ~ 300 GeV for all values of tan /3 (barring unnatural cancellations in loops between the charged 
Higgs and new light superpartners) |S2]. In the L2HDM the coupling to all quarks can be made arbitrarily small as 
tan/3 is increased. For large tan/3, the loop-induced contributions of to tau lepton decay become important [63) . 
but as shown in Figure]^ the available parameter space for charged Higgs in the L2HDM is significant. Figure 2 in 
Ref. [62] shows the analogous plot for the Type H 2HDM. Comparing the two it is evident that the L2HDM allows 
for a much larger parameter space than the Type H case. 

As the charged Higgs couples to the Standard Model-like Higgs, it will contribute to the loop-induced /177 decay. 
This coupling is proportional to A 3 , and so for arbitrary L2HDM models can significantly decrease the diphoton width 
of the 125 GeV Higgs |48|. However, for MSSM-like couplings as in Eq. (§, this loop contribution is small, resulting 
in only a ~ 3% decrease in the diphoton decay width for a charged Higgs with mass at the LEP-H limit. This is in 
mild tension with the measured value, as the CMS and ATLAS measured h — > 77 rate is higher than the Standard 
Model value unmi]. 


C. Muon g-2 

As seen in the limits on the charged Higgs, the presence of light particles in the spectrum can contribute to 
Standard Model processes at the loop-level. While in almost all cases, the spectacular success of the Standard Model 
in matching observations means that these loop contributions can only be used to set limits on new physics, the 
measured value of {g — 2)^ differs from the Standard Model prediction by over 3cr: Aa^ = (262 ± 85) x 10“^^ [iTZ] 
where = {g — 2)^ /2. New pseudoscalars and scalars introduce contributions to {g — 2)^ at both one- and two-loops, 
which might help account for this discrepancy. We calculate these corrections as a function of mediator mass and 
tan/3 following Eqs. (8 — II) in Ref. |S7|, for both a CP-even scalar H and CP-odd pseudoscalar a in the L2HDM. 

The loop corrections from the scalar or pseudoscalar are proportional to its coupling with the muon. Since the 
pseudoscalar a couples to leptons proportional to tan /3 in the L2HDM, the contribution from pseudoscalar can greatly 
enhance {g — 2)^. In Figure Isl we plot the contribution of the scalars and pseudoscalars to as a function of scalar 
mass for various values of tan/3. Pseudoscalars contribute positively to the deviation from the Standard Model value, 
while scalars contribute negatively for large tan /3. We see that a light pseudoscalar with high tan /3 can account for 
the anomalous measurement. Note that the scalar contribution is negative, however at high masses this contribution 
is negligible. However, for the L2HDM to be successful in explaining the observed {g — 2)^ anomaly, a large mass 
splitting between the light pseudoscalar and the scalars is necessary. Though this is possible in a general L2HDM, 
it is difficult to accomplish in supersymmetric models without large loop corrections. Our results agree with those 
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FIG. 2: Constraints on charged Higgs mass mH± in the LH2DM from LEP II results [SS], r decay |B3], top decay 
[Sn], and B —>• (calculated using the procedure described in Ref. 


found in the recent work Ref. [IS] . 



FIG. 3: Deviation of from Standard Model value (Ao^) due to L2HDM scalars and pseudoscalars as a function of 
scalar mass for various values of tan/3. Also shown is the experimental best-fit region m- 
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D. Direct Collider Searches 

As with any proposed new particle with electroweak-scale masses, searches for heavy neutral bosons in extended 
Higgs sectors have been a focus at colliders. Unfortunately, since we now know we live in the limit sin(/3 — a) « 1, 
many processes which would contribute significantly to production of the additional CP-even scalar H in a generic 
2HDM are now known to be highly suppressed. This is because the couplings to gauge bosons of the H is proportional 
to cos(/3 — a), so the WjZ + H associated production and vector boson fusion mechanisms are effectively eliminated 
in the alignment limit. 

The production of the CP-even H in association with the pseudoscalar a through the ZHa coupling is not suppressed 
in the alignment limit, and is in fact maximal. Searches for neutral Higgses done by LEP-H looked for neutral boson 
production from the process e+e“ —>■ Z* —?> H{h) a, assuming each decays into pairs of tans (as will be the dominate 
decay in the L2HDM). This search concluded that there cannot exist both a light pseudoscalar and a light scalar 
which are kinematically accessible at LEP-H [55]. Roughly speaking, the sum of H and a masses are required to be 
greater than 180 — 200 GeV. At this time, this LEP-H bound places the only significant limit on the mass of the H 
and a, given that the L2HDM must reside in the alignment limit. As with resolving the {g — 2)^ anomaly, we find 
that very light scalars or pseudoscalars are possible, but not both. 

The production of the H or a via gluon-fusion is also suppressed in the L2HDM (assuming tan (5 > 1), as this vertex 
is induced via top-loops and is therefore proportional to cot fi. Associated production via top or bottom quark pairs 
is similarly reduced, as was the case with as discussed previously. One might consider associated production of H 
and/or a alongside tau pairs, but large production cross sections can only be obtained for very light nia or mu, and 
for tan/3 so large that tan/3 x rrir/v > 1. This leaves the associated production of H and a through the unsuppressed 
coupling to Z as the main production mechanism at the LHC. For tan/3 > 3, the decay to tau pairs dominates over 
other fermionic channels for both the H and the a, so both of these production mechanisms would result in four tau 
leptons in the final state. If ttIq < TOij/2, the channel iL —>■ aa can have a large branching ratio, resulting in six-tau 
final states. 

The multi-tau final state falls under the multilepton search category at the LHC |69|. We simulate production 
of L2HDM scalars and pseudoscalars using MadGraphS [751 ITT] and PythiaG, and DelphesS m for detector 
simulation. We find that the existing CMS multilepton search performed on 20 fb“^ of 8 TeV data does not yet have 
sufficient sensitivity to place limits on the production of leptophilic Higgs scalar/pseudoscalar pairs with masses at 
the LEP-H threshold. The 13 TeV Run-H will have signal production cross sections larger by a factor of two, but 
with similar increases in the primary irreducible dibosons backgrounds. Therefore, discovery of a leptophilic H and 
a at the LHC would require additional improvements in the multi-tau final states in order to extract the small signal 
from the background (see e.g. Ref. [75]). 


III. DARK MATTER 


We now consider the addition of dark matter to the leptophilic models. Both scalar [74| [^ and fermionic m 
dark matter have been considered in this context previously, but as in the broader consideration of the L2HDM, the 
discovery and measurement of the 125 GeV Higgs suggests a reassessment of the possibilities with the new evidence in 
mind. In this work we assume fermionic dark matter, as we are motivated by possible connections of L2HDM models 
with supersymmetry, where the dark matter would be fermionic neutralinos. In our non-supersymmetric L2HDM, we 
must add dark matter in by hand, by introducing two SU{2)i^ doublets with hypercharge -1-1/2 and —1/2, as well as 
a gauge singlet (these fields are analogous to higgsinos and binos in a supersymmetric model): 


•01 = 




V's- 


Charging the 0i as odd under the same Z 2 as in the Higgs potential, and allowing only soft-breaking terms, the dark 
matter-Higgs Lagrangian for the neutral states then takes the form 

-Cv- 2 -2/l'010L0S - y2020Q0S - ^msi’si’S - 1*0102 + h.C. (11) 

The couplings yi and 1/2 are, in a general theory, free parameters. For specificity, and in order to draw a connection 
with supersymmetric models, we will assume yi = — 1/2 = g' ^ 0.35 for numeric calculations, as would be the case for 
supersymmetric higgsinos and binos coupling to the Higgs sector through the superpotential. A full supersymmetric 
version would give rise to higgsino and gaugino-matter couplings with Yukawa couplings oc g'. After EWSB, this gives 


a mass matrix for the neutral components of 
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( 12 ) 


This matrix can be diagonalized by a unitary transformation. We relate the original fields to the mass eigenstates 
with the coefficient matrix: 


/ Nil Ni2 Ni3 \ / Xi \ ( 

N21 N22 N23 ^2 = ^2 

\ -^31 N 32 N 33 ) \xz ) \'<Ps J 


where the eigenstate xi = X corresponds to the lightest neutral particle, and thus is a dark matter candidate. Note 
that the coefficient /r will set the mass of the charged states. To satisfy LEP-II constraints on new charged particles, 
we take p, > 100 GeV for the remainder of the paper. 

In Dirac notation, the dark matter x has interactions with the physical Higgs fields and the Z boson of 


N 31 _ N 31 

^int D Hxx-^ {CaVlNii + Say2N2l) - hxx-^ {SaVlNii - Cay2N2l) 


(13) 


{sfsyiNii + cpy2N2i) - [liVnp - jiVzip] {xi^l^x)Zt,. 
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These terms define the dark matter couplings ghxx^ 9 Hxx^ 9axx ^'^d gzxx- that in the alignment and large 

tan/? limit, both the H and a couple primarily to the Nn component (the “leptophilic higgsino,” if viewed as a 
supersymmetric model), while the h couples mostly to the N 12 component (the “quark-like higgsino”). 

As the dark matter arises from mixing SU{2)l doublets with singlets, it couples to both the 125 GeV Higgs and 
the Z, and so can contribute to the invisible width of these particles. The direct decays have widths 


r(/i ^ xx) 


r( 2 ’ ^ xx) 


Stt y 



3/2 


5 I 


xx^^ 

2Att 



(14) 

(15) 


As discussed previously, the upper limit of the invisible branching ratio is 54% of the total Higgs width [5D1154j . 
The upper limit on the total Higgs width is ^ 17 MeV [50], and so at most we must require T/n^^^ < 9 MeV. Assuming 
that, in addition to the dark matter interaction, all the h couplings are exactly those of the Standard Model, then 
the limit on the invisible branching ratio translates to T^.^^ < 5 MeV. The uncertainty in the Z width constrains 
^ 2 MeV [77]. The invisible widths for Z and /i as a function of dark matter mass are shown in Figures]^ 
and [^ respectively, assuming fi = 200 GeV. As increases, the SU{2)l component of the dark matter decreases, 
and so the constraints are more easily satisfied. Note that rriy^ cannot be made arbitrarily small for fixed values of 
the Yukawas yi,y 2 and /x, thus explaining why some lines in Figures [^ and [^ are cut off at the low mass end. For 
our choices of parameters shown in these figures, the ghxx coupling will receive a cancellation from the Y = -1-1/2 
and Y = —1/2 components of the dark matter. This cancellation occurs at lower for higher tan/3, explaining the 
cancellation in the tan/3 = 100 curve in Figure]^ 

In addition to direct decays of h to dark matter, h could decay to a pair of pseudoscalars, which themselves can 
decay to dark matter pairs (assuming that each decay is kinematically allowed). The constraints on the h ^ aa 
discussed in the previous section apply here. One might also consider whether invisible decays of the a and H to dark 
matter could allow the LEP-H bounds on the production e~e~^ —>■ Z* —>■ Ha [5S] to be evaded, by turning a four-tau 
final state into one of missing energy only. However, this is possible only for very low values of tan /3, as the branching 
ratio of a and H to tau leptons increases at large tan/3, while the coupling to dark matter asymptotes to a constant 
in this limit. In Figure]^ we show the upper limit on tan/3 from LEP-H as a function of scalar and pseudoscalar 
mass, for 1 GeV dark matter with the maximal doublet component which is in agreement with the invisible width 
constraints. This choice of parameters maximizes the branching ratio of the a and H into dark matter and so most 
efficiently avoid the constraints from the tau search at LEP-H. 
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FIG. 4: Invisible Z width as a function of dark matter mass for various values of tan/3 with ^ = 200 GeV. The limit 
is ^ 2 MeV [77]. 



FIG. 5: Invisible Higgs width as a function of dark matter mass or various values of tan f3 with /i = 200 GeV. The 
limit derived from the maximum Higgs width is Th^xx ^ ^ MeV |50|, while the limit derived from the SM Higgs 
width is < 4.8 MeV [SOHM] . 


A. Direct Detection 

Direct detection experiments measure the recoil energy from WIMP-nucleus scattering, placing an upper limit on 
the dark matter-nucleon elastic scattering cross section. This allows us to place bounds on the parameters that make 
up the dark matter-scalar and quark-scalar couplings. Dark matter-nucleus scattering is dominated by mediation from 
the GP even scalars. The pseudoscalar mediated interaction is velocity suppressed and so has negligible contribution 
to the scattering cross section. Dark matter-nuclei interactions mediated by scalars induce a spin-independent cross 
section which is constrained by a number of experiments. Presently, the strongest bounds come from the LUX 
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(GeV) 

FIG. 6 : Upper limits on tan jS from LEP-II constraint [ 68 ] assuming decay into either dark matter or rr, setting the 
mass of dark matter to 1 GeV and the doublet component Nn = 0.2. 


experiment [75] for > 6 GeV and by GDMS-lite at lower dark matter masses m- 

The fundamental Lagrangian parameters are translated into dark matter-nucleon scattering cross sections using 

m 


^v—-D.n — 


^ .2 
_ J p.n 5 
TT ^ 


I 2 p^n 1Tlp,n ^ 

« ma 27 ^™ 2 ^ 


fp,n — ^ 

q=u,d,s ^ q—c,b,t 

Sg — jy^2 9Hqq dHxX + J^2 9hqq 9hxx 


m„ 


(16) 

(17) 

(18) 


where /x is the dark matter-nucleon reduced mass, is the effective dark matter-quark coupling, and the parameters 
/P’” and /^Q are proportional to the quark expectation operators in the nucleon. These must be extracted from lattice 
QCD simulations [5TI - I5H] . and we adopt the values from Ref. |5S]. For the purposes of this paper there is no signihcant 
difference between the proton and neutron fp^m and so our dark matter scattering is essentially isospin-conserving. 

In Figure we show the direct detection cross for various p, and tan (3 values as a function of dark matter mass, 
with iTiH = 500 GeV (though note that, due to the leptophilic nature of this scalar, a lower mass for H does not 
significantly change the result). We see that there is space to accommodate LUX bounds, especially as increases. 
Note again two accidental cancellations which can reduce the measured cross section. One decreases of the SU (2) l 
doublet components’ contributions to the cross section, which shifts to lower dark matter masses as tan jS is increased. 
The second occurs as —>■ fi resulting in 7 Vi 3 —)■ 0, which we have chosen to be 200 GeV here. 


B. Indirect Detection 

The L2HDM with fermionic dark matter model allows for the annihilation of dark matter into Standard Model 
particles in the Universe today. Indirect detection experiments search for unexplained excesses of gamma rays, 
positrons, or other Standard Model final states coming from areas of high dark matter density. Whereas direct 
detection experiments are sensitive to s-wave scattering mediated by scalars, thermally averaged annihilation cross 
sections mediated by scalars are p-wave suppressed, i.e. proportional to . The velocity of dark matter today is very 
small (< 10“^c), and so p-wave annihilation is not detectable. However, annihilation through a pseudoscalar has a 
velocity independent component and therefore is the main contribution to {av). The thermally averaged cross section 
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FIG. 7: Spin-independent cross section as a function of dark matter mass for various values of /r and tan/3. Limits 
from LUX are given in green. We take mn = 500 GeV. 


for {xx //) takes the form: 


{av) 


X"/ 

/ 


[{m^ — 4 to ^ r+min] 



3m} 



(19) 


where Nf is the number of colors of the fermion /, and T is the temperature of the dark matter. We can safely 
take T = 0 since the temperature today is very small. As the a couples primarily to tau leptons in the large tan /3 
limit, which decay into final states with large number of photons, the primary indirect detection search channel will 
be gamma rays. 

Areas of expected high density dark matter in our Galaxy can be found in nearby satellite galaxies [HI [5M5H] , 
extra-Galactic sources [551 - 155] . or the Milky Way Galactic Center. A notable excess of gamma rays coming from 
the Galactic Center has been reported by various analyses of Fermi Gamma-Ray Space Telescope (FGST) data [43l - 
146] . The source of these gamma-rays is uncertain but can be interpreted as coming from dark matter annihilation. 
Annihilation to taus is one of the final states which fit the excess spectrum reasonably well [5111105] . when the dark 
matter mass is in the range ~ 9 — 11 GeV. In Figure [^we show the thermally cross sections as a function of tan /3 for 
various values of /i and ma, fixing = 10 GeV. The shaded regions are possible values for {av) in our model. As 
expected the cross section decreases as the mass of the pseudoscalar increases, however we can see that there is still a 
large range of parameters which fit the Galactic Center excess, even when the pseudoscalar and dark matter masses 
are far from being in resonance. It is notable that a leptophilic Higgs model furnishes a viable mediator between 
the visible sector and the dark matter which provides a sufficiently large annihilation cross section; in Type-II Higgs 
models, a Higgs mediator which can explain the Galactic Center anomaly without a resonance enhancement tends 
to require additional light states which are well-constrained from direct detection and collider searches |106H112] . As 
the mediators considered here are leptophilic, such constraints are easily avoided. 


C. Thermal Relic Abundance 

Measurements of CMB anisotropies from telescopes such as Planck give the dark matter component of the Universe’s 
energy density to be = 0.1187 ± 0.0017 [113j . The standard Boltzmann relic density calculations imply the 
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tan/3 


FIG. 

GeV 

from 


8 : Thermally averaged annihilation cross section as a function of tan /3 for ma = 25 GeV (blue) and nia = 50 
(green), assuming = 10 GeV. Also shown are best fits to the Galactic Center excess in the rr final state, 
Refs. [Tomrm5] . 


corresponding thermal annihilation cross section for dark matter of {av) ~ 3 x 10“^® cm^/s. If we assume that the 
dark matter is a thermal relic we can calculate the necessary couplings to give the correct thermal annihilation cross 
section. The calculation is the same as in the case of indirect detection except that we evaluate T at the freeze-out 
temperature Tf = [80]. 

Along with showing available parameter space for indirect detection, Figure [^ also demonstrates the region of 
parameter space with the correct thermal abundance obtained through couplings with the pseudoscalar. We see 
that the canonical value of the thermal annihilation cross section can be achieved for many values of tan/3 and rua, 
although high values of ma require increasingly higher values of tan /3 as expected. 

The best fits for the Galactic Genter anomaly do not coincide with the thermal relic cross section, a concern that our 
model shares with other explanations of the excess assuming annihilation to taus. It is possible that additional p-wave 
processes contributed to the thermal relic cross section which are simply inactive today due to velocity suppression. 
Alternatively, processes going to final states which do not contribute significant gamma-ray flux in the Galactic Center 
(for example, annihilation to neutrinos) can serve to boost the thermal relic cross section for freeze-out in the early 
Universe. Our results shown in Figure should be taken to demonstrate that thermal cross sections of the correct 
magnitude for both thermal relics and the Galactic Center can be relatively easily obtained in the L2HDM with 
fermionic dark matter, without significant fine-tuning. 


IV. LEPTOPHILIC SUPERSYMMETRY 

Throughout this paper we have often set unknown couplings in the Higgs and dark matter sectors to their “supersym¬ 
metric” equivalent in order to both reduce parameter space and to make contact with the most familiar interpretations 
of two-Higgs doublet models, as realized in the MSSM. However, the L2HDM cannot be supersymmetrized without 
adding new Higgs fields, as the Yukawa interactions for the Standard Model fermions must come from a holomorphic 
superpotential, forbidding down quark masses from the d>Q term in Eq. 0. The minimum number of Higgs doublets 
required for a leptophilic supersymmetric model which contains the Standard Model is four, i7„, and Hn, 
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with the down- and lepton-Higgses having hypercharge —1/2 and the up-type and Hn having -1-1/2. This model was 
introduced in Ref. and we adopt much of their notation. 

Of the four Higgs doublets, we assume as in the L2HDM that the and i?„ are odd under some Z 2 , along with 
the lepton doublets. Then the superpotential for the theory is (again suppressing flavor indices) 


W = VuQlUbHu + VdQLdnHa + ygLLenHg + 


( 20 ) 


The soft supersymmetric-breaking potential will break the Z 2 symmetry, 

Hsoft = AW + A\Hd? + + AnW + {AlHuHd + A 2 Hnm + AHuHi + AiHnHd + h.C.) . ( 21 ) 

The quartic couplings are provided by supersymmetric H-terms, as in the MSSM. We assume that each Higgs 
obtains a CP- and 17(l)em-conserving vev, with A- analogy with the ratio tan/3, we define 

three angles: 

tan^ /3 = , tan 7 = —, and tan 7 ^ = —. (22) 

vj + Vd Vl 


As in the L2HDM, the angle (3 describes the amount of the EWSB vev which lies in the lepton sector, and we are 
primarily interested in the large tan/3 limit. We refer readers to Ref. |114j for the mass matrices of the Higgs sector 
after EWSB. 

As in the MSSM, the four-Higgs supersymmetric model has a natural dark matter candidate, the lightest neutralino, 
here made up of winos W, binos B, and neutral higgsinos: hd, hn and hi. The mass matrix for the dark matter 

(in the hu, hd, hn, hi, W, B basis) is 


M^ = 


( 0 

Ml 

0 

0 

gvs 

", 

\ 


Ml 

0 

0 

0 

gvspc^^ 

,V 2 


0 

0 

0 

M2 

gvc^ 


9 g vc(,s^^ 


y/2 


V2 


0 

0 

M2 

0 

gvc0 

,C 2 

g vcqc~^^ 

72 


gvs0S~^^ 

g vs0S~y^ 

V2 

±2 

gvsi 3 C~y^ 

g vsi^c~y^ 

V2 

V2 


g vc0S~y^ 

±2 

,±2 

gvc^ 

g vcjsCy^ 

\/2 

V2 

M 2 

0 

0 

Ml 


(23) 


/ 


This matrix is diagonalized by the 6 x 6 matrix N. Assuming Mi ^ yi, fi 2 ^ M 2 , the lightest dark matter candidate 
will be primarily B with an admixture of higgsino. The couplings of the neutralinos to the Z are: 


^ _ |/V,,d|2 - \nW WAWi.- 


(24) 


The couplings to the W CP-even and CP-odd Higgses are: 


'^int — 2 [9^k,W 9 ^kjS] Ui^d^k^d T Ui^n^k,n 


(25) 


+ 2 — gNk.s] [{Vt^uNk^u — Vi^dNk.d + Vi^nNk.n — Vi/Ni^l)] hiXkXk, 


where U and V are the rotation matrices for the CP-even and CP-odd Higgses hi and Oi. 

After EWSB, the physical scalar sector has four CP-even scalars hi, three CP-odd pseudoscalars Ui, and three 
charged scalars hA As the Higgs quartics are set by the gauge-interactions of SU(2)l doublets with hypercharge 
±1/2, the CP-even Higgs scalar that has Standard Model-like couplings to the W and Z bosons has a tree-level mass 
bounded from above by mz- Therefore, large loop-contributions to the mass of the Standard Model-like Higgs are 
still required, as in the MSSM, and the Little Hierarchy Problem [15] remains. 

However, as in the L2HDM, this four-Higgs model does allow for the possibility that the 125 GeV Higgs is not the 
lightest Higgs scalar, unlike in the MSSM. As discussed in Section [Hj3, in a Type-H 2HDM, as in the MSSM, the 
charged Higgs cannot be lighter than ~ 300 GeV without requiring light superpartners who’s contribution to 6 —> 37 
is fine-tuned to cancel the charged Higgs loop. In supersymmetric models, the mass of the charged Higgs, GP-odd 
pseudoscalar, and one GP-even scalar are set by a common mass scale with relatively small splittings. Thus, the high 
mass of the charged Higgs requires the Higgs at 125 GeV to be the lightest in the MSSM. 
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Given the large number of possible parameters and the need for large loop corrections to obtain a 125 GeV aligned 
Higgs (as in the MSSM), we do not perform scans to ht to data, but restrict ourselves to a general discussion of the 
possibilities. In the large tan /3 limit of the four-Higgs supersymmetric model, the mass matrices for the Higgs sector 
become nearly block diagonal (barring very large ^3 and ^4 terms). As a result, the GP-even scalar with Standard 
Model-like couplings to the W and Z bosons (which we must identify with the 125 GeV Higgs discovered at the 
LHC) is primarily composed of an admixture of and Hd, as are the pseudoscalar and charged fields eaten by 
the Z and W bosons. The remaining and Hd states, consisting of a CP-even neutral scalar, a CP-odd neutral 
pseudoscalar, and a charged scalar, must be heavy enough to avoid large flavor-changing decays in the quark sector, 
as in the MSSM. The leptophilic components, two each of the neutral scalars, pseudoscalars and charged pairs, can be 
a generic admixture of Hi and iL„ fields, though if tany^ ^ 1 or ^ 1 , this will tend to result in each mass eigenstate 
being a nearly pure Hi or state. At tree level, two triads composed of a scalar, pseudoscalar, and charged Higgs 
will cluster in mass, with the scalar and pseudoscalar masses very close together, and the charged Higgs heavier by a 
addition to its mass squared. 

These leptophilic Higgses can themselves be very light. The charged Higgses composed of H^ and H^ must be 
heavier than the LEP-H limit of ~ 92 GeV, but as seen in the L2HDM have very few other meaningful constraints. 
The LEP limit can be easily satisfied, even if the neutral components of the triad is well below mz- 

The GP-even and -odd components of a triad can be produced through their coupling to the Z, as in the L2HDM. 
Therefore, the LEP-H limits on e“e+ —)■ Z* —)• Ha — )■ 4t, as discussed in Section |Hj3 apply here. As before, the sum 
of the scalar and pseudoscalar masses must be > 180 GeV. Barring large loop corrections splitting the scalar and 
pseudoscalar, this effectively places a lower limit on the mass of these new Higgs pairs of ~ 90 GeV. 

However, given that large loop corrections must be applied to the Standard Model-like Higgs to lift its mass to 
125 GeV, it is not implausible that the loop corrections from the matter sector can introduce large additional quartics 
to the leptophilic sector as well, splitting the mass of the scalar and pseudoscalar components. It is interesting to note 
that, in the L2HDM with supersymmetric-like couplings, we found that the pseudoscalar could not be below half of 
the 125 GeV Higgs mass without introducing non-standard decays that are experimentally ruled out. This situation 
can be mitigated in the supersymmetric case if the lightest pseudoscalar is an approximately equal admixture of a„ 
and ai- In this case the coupling to the Standard Model-like Higgs is suppressed due to a cancellation between the 
two components with opposite hypercharge. This does allow a leptophilic pseudoscalar in the 60 GeV mass range 
preferred by {g — 2)^ and suggested by the Galactic Genter anomaly, though we again stress that this requires a large 
splitting in the scalar-pseudoscalar masses, introduced by loops. 

The results from the L2HDM with fermionic dark matter added by hand can be ported to the supersymmetric 
version, where the dark matter sector arises naturally. A primarily-bino dark matter particle with 0(10%) admixture 
of higgsino components can easily be obtained when Mi jli M 2 . Large tan/3, combined with a light pseudoscalar 
from the Hi and H^ doublets will result in annihilation cross sections consistent with the Galactic Center anomalies or 
thermal annihilation, while the direct detection constraints are easily satisfied due to the decoupling of the quark-like 
Higgses. As in the L2HDM, the light scalars and pseudoscalar mediators between the dark matter sector and the 
visible sector have extremely weak constraints from collider searches; as a result, the most easily obtained signature 
of this leptophilic supersymmetric model would be through the dark sector. 


V. CONCLUSIONS 

Post-Higgs discovery, one of the high priority tasks at the LHC is to determine the full structure of the Higgs 
sector which is responsible for electroweak symmetry breaking. While the Higgs discovered at 125 GeV appears very 
consistent with the predictions of a single Higgs doublet in the Standard Model, our measurements are still only 
accurate at the ^ 10 % level for couplings to the electroweak gauge bosons, and significantly weaker when it comes to 
direct measurements of the couplings to fermions. Much work remains to be done. 

Two-Higgs-doublet models are a simple and attractive extension to the Standard Model Higgs sector. However, 
much of the experimental and theoretical effort has been dedicated to a particular version of such models, the Type-H 
model as found in the MSSM. Strong indirect limits can be set on such models, due to the coupling of both Higgs 
doublets to the quark sector. In particular, the charged Higgs must be heavy, which in the MSSM indicates that the 
discovered Higgs is the lightest of the physical scalars and the remaining Higgs are significantly decoupled. 

In the light of the LHC results, we reconsider an alternative two-Higgs model, in which one doublet couples to both 
up- and down-type quarks while the other couples to the leptons only. This model easily evades most experimental 
searches for additional components of the Higgs sector, especially in the large tan /3 regime where the new mass 
eigenstates couple mostly to leptons, rather than quarks. Additional Higgs scalars can be lighter than the 125 GeV 
CP-even Higgs, with the strongest constraints coming from LEP-H. These limit the charged Higgs to be heavier than 
^ 92 GeV, and require the sum of the additional CP-even and CP-odd Higgs masses to be greater than 180 — 200 GeV. 
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Given the relative alignment of the 125 GeV Higgs, should leptophilic Higgses exist, their discovery at the LHC 
will be difficult. From the LHC Higgs measurements, we now know that the most promising channels, considered 
previously munn] would have small cross sections. Direct production of the scalar and pseudoscalar particles 
through VBF or gluon fusion process is suppressed by cos(/3 — a)^ or cot^ /3. The pair production mechanism of a 
scalar and a pseudoscalar through the Z followed by decays to tau lepton pairs is not suppressed, but is an extremely 
difficult search given the efficiency for tau-tagging and the large backgrounds. For example, we have demonstrated 
that the existing multilepton analysis from CMS in Run-I [5S] is not sensitive to scalar/pseudoscalar pairs with masses 
at the LEP-H bound. The increased cross section at Run-H will be helpful, but a full analysis of the tau backgrounds 
will be necessary and the channel will remain difficult. New ideas may be necessary. 

Indeed, the discovery of leptophilic Higgses might be most easily achieved through its coupling to dark matter. 
Motivated by supersymmetric extensions of the leptophilic model, we consider fermionic dark matter composed of 
SU{2)L doublets mixed with singlets - analogs of higgsino-bino dark matter. In the large tan /3 limit, the dark matter 
has significant interactions with the leptophilic Higgs scalars, through the doublet component. We hnd significant 
parameter space for thermal dark matter, with possible masses extending below half the mass of the 125 GeV Higgs or 
the Z mass without violating the invisible width constraints on those particles. Dark matter indirect detection rates 
are naturally suppressed compared to the expectations from the Higgs portal in Type-H model. However, barring 
accidental cancellations, the direct detection rate is in reach of the next generation of direct detection experiments. 
The large couplings to the leptophilic pseudoscalar allow for significant s-wave annihilation in the Universe today, 
allowing for a 10 GeV dark matter particle annihilating into tau leptons, as has been suggested as a possible fit to 
the Galactic Genter anomaly with mediator masses < 90 GeV. 

We revisit the supersymmetric extension of this model, including four Higgs doublets, which is leptophilic in the 
limit where the majority of the EWSB vev resided in the quark sector. As in the MSSM, significant loop-corrections are 
required in order to bring the Standard Model-like Higgs up to the observed mass. Thus, though the supersymmetric 
leptophilic model allows the interesting possibility of new Higgses lighter than the one already discovered, it does not 
address the Little Hierarchy Problem. The attractive properties of dark matter in a leptophilic two-Higgs model are 
found as well in the supersymmetric version: low mass mediators capable of fitting the Galactic Genter anomaly would 
not have been seen in colliders, and unlike in the MSSM, would not induce large direct detection signals in tension 
with experimental results. Mixing between the two leptophilic Higgs doublets can also allow a pseudoscalar lighter 
than half the 125 GeV Higgs mass, which is not possible in the leptophilic two-Higgs model. As in the L2HDM, the 
discovery of the extended Higgs sector by itself is difficult at the LHC, as the production cross sections most conducive 
to LHC searches are suppressed in a leptophilic Higgs model with one nearly-Standard Model-like Higgs, as is the 
experimental situation we find ourselves. The addition of supersymmetric particles allows for larger production cross 
sections, as colored particles will undergo cascade decays into the neutralino/chargino sector. These cascades would 
preferentially include tau leptons in the supersymmetric leptophilic model. Though searches including tau final states 
have been performed, they are experimentally challenging, and existing constraints are weak and the potential for 
improvement is large. 
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